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We consider the application of endpoint techniques to the problem of mass determination 

n 

for new particles produced at a hadron collider, where these particles decay to an invisible 
qj particle of unknown mass and one or more visible particles of known mass. We also consider 



decays of these types for pair-produced particles and in each case consider situations both 
with and without initial state radiation. We prove that, in most (but not all) cases, the end- 
point of an appropriate transverse mass observable, considered as a function of the unknown 
mass of the invisible particle, has a kink at the true value of the invisible particle mass. The 
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co-ordinates of the kink yield the masses of the decaying particle and the invisible particle 



We discuss the prospects for implementing this method at the LHC. 
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^ I. INTRODUCTION 

With the imminent start-up of the Large Hadron Collider (LHC), and the hope that this machine 
could produce the first man-made dark matter, the question of how to measure the masses of heavy, 
invisible particles at hadron colliders has never held more importance. 

Cho et al. [I] have recently made a remarkable claim regarding a method for measuring invisible 
particle masses. They studied identical pairs of three-body superpartner decays (g — > qqxi)i with 
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one of the daughter particles being invisible, in the special case where the sum of the transverse 



momenta of the parents vanishes. They identified "kinks" (see, for example, Fig. 10) in the graphs 
of suitably defined transverse- mass variables, considered as a function of the (a priori unknown) 
mass of the invisible daughter. The importance of their claim is that it provides, in principle, 
a method to determine both of the unknown masses in the problem - namely the mass, mo, of 
the decaying parent particle and the mass, mi, of the invisible daughter. These masses are given 
simply by the co-ordinates, (mo, mi), of the kink. 

In [2] , one of us gave a proof that a similar kink occurs even in events of a much simpler type, that 
is, in events in which a single parent with a non-zero transverse momentum distribution undergoes 
a two-daughter decay. It was claimed in j2] that such kinks should be a generic phenomenon. 

We need hardly stress how useful an "transverse mass kink method" for determining sparticle 
masses would be, if experimentally viable. In particular, if the invisible particle were indeed 
the dark matter particle, then knowledge of its mass would be a boon for the astrophysics and 
cosmology communities. What is more, knowledge of the spectrum of, say, superpartner masses, 
would be of much help to those attempting to discover the method by which supersymmetry 
breaking is mediated to the Standard Model. 

When the results of [l j and |2j first became public, it was not yet clear that there was a 
difference in the manner in which the kinks arise. At first sight, both papers reported kinks in 
related transverse variables. But, as we go on to show in this paper, the sources of the kinks in [1] 
and [2] are not the same. The results have no direct relevance to each other, although each paper 
separately provided very important results. 

One of the key purposes of this paper is to explain the connections between the results of [1] and 
j2] and the relevance for the LHC. While Cho et al. dealt with the class of events where the sum of 
the two parents' transverse momenta is strictly vanishing, [21 dealt with the class of events in which 
the transverse momenta could be arbitrarily large. Neither of these situations is a particularly good 
approximation at a hadron collider, since a realistic sample of events will be somewhere between the 
two - they will have a distribution of possibly small, but certainly non-zero, transverse momenta. 
We address this more physically relevant case and find that a reasonably sized sample of realistic 
events will still produce a kink at (mo, mi), from which, in most cases, we can hope to extract the 
masses. 

We discuss the implications these results have for experimental observability and explain how 
the events containing (pairs of) three-body decays considered by [1] were able to generate kinks 
even without approaching the asymptotic extrema examined in [2 . We comment on how the kinks 



seen in [T] would be modified by a (physically reasonable) admixture of events in which the parents 
have significant transverse recoil. We also point out that if nature were to provide us with enough 
events far above threshold, then pairs of two-daughter decays would generate a kink, even though 
the mechanism by which the kink is generated is independent of the arguments in [1] . 

In doing so we extend the analyses of both [I] and [2] . We analyze the general case of N + 1- 
daughter decays, including cascade decays. We also examine cases with identical pairs of such 
decays/cascades. We consider events both with and without recoil of the parent particle(s) against 
a secondary system of significant transverse momentum (respectively SPT and ZPT). We also 
perform Monte Carlo simulations of a variety of different examples. Since the list of cases studied 
is not short, the reader may find it useful to refer to Fig. [T] and Fig. [2] which show and label various 
scenarios diagrammatically. 

The layout of the paper is as follows. In Section O we (re)introduce endpoint analyses as 
methods of measuring masses. We discuss the extensions of the transverse mass variable, tut, to 
cases where parent particles are pair-produced, and where the mass of the invisible daughters is 



unknown. In Sections III A and III B we consider single parent decays (both SPT and ZPT) and 
determine the complete set of extrema of an appropriately-defined my. We analyze the nature of 
the various extrema, identify the global maximum, and examine the fraction of events which might 



be expected to lie near that maximum in each case. In Section III C we consider events containing 



pairs of identical decays of the type discussed in Section III A In Section IV , we present the results 
of some Monte Carlo simulations which illustrate these results. We conclude in Section [V] with 
a summary of our findings, a comment on the relevance of SPT versus ZPT for the LHC, and a 
discussion on the measurability of such edges in LHC-like scenarios. 

II. ENDPOINT TECHNIQUES 

The determination of mass in the complex environment of hadron collisions is not usually easy. 
For example, backgrounds from known physics are rather large, misidentification of particles and 
jets is a problem, and we have very little direct information concerning the longitudinal velocity of 
the centre-of-mass frame of the two primary interacting partons. 

These difficulties are compounded if the new particles are such that they decay with significant 
branching fraction to particles that interact sufficiently weakly as to be invisible as far as the 
detector is concerned. In any one event, these invisible particles carry off kinematic information, 
namely their momenta, and also their energies, if their masses are unknown. 



The loss of these particles means that we cannot always measure particle masses on an event- 
by-event basis. But since we are able, in principle, to observe multiple events, the situation is not 
beyond hope. The historical example of the discovery, and subsequent mass determination, of the 
W-boson, is rather instructive in this respect. 

The VK-boson was first discovered in the UA1 and UA2 experiments [31111 [5], through its sizable 
leptonic decay to an electron and a neutrino. The latter is, of course, invisible in the detector, and 
kinematic information thus is lost in each decay. Nevertheless, the transverse mass, defined by 

m T = m l + m l + 2(e e ej / - p e .Pi,), (1) 

where e e = ym^+p*, is observable. The electron mass, m e is known a priori, and the neutrino 
mass, m u , is negligible. The transverse momentum of the electron, p e (and its magnitude, p e ), can 
be measured directly, and the transverse momentum of the neutrino, p y , can be inferred from the 
missing transverse momentum in the event. Now, for a single event, m^ is not of great interest; its 
importance lies in the fact that, in the limit of a narrow width for the W-boson, m^ is bounded 
above by the mass of the W, mw, and is equal to mw m events where the electron and the neutrino 
have the same rapidity. To see this, note that the four-momenta, (E,p,q), obey the constraint 

rriiY = ml + ml + 2(E e E u -p e .p u - q e q u ). (2) 

Now, the rapidity is given by 

^^E—q' (3) 

and the vanishing of the relative rapidity r] e — r\ v implies 

E e q u = E u q e , (4) 

whence 

E e E u - q e q u = e e e u , (5) 

and 

m\v = mr- (6) 

Thus, by computing tut for many events, one can obtain a distribution of ttt-t values whose upper 
endpoint is close to the true value mw- Indeed, fitting the tut distribution using data from CDF 
provides the most precise single direct measurement of the W mass [6] . 



Similar and more general endpoint techniques are expected to be useful at the LHC. Extremal 
values of the kinematic observables typically contain information about the particle masses, as 
these are the parameters which determine shape of the boundary of multi-particle phase space. 

Although endpoint techniques suffer from the disadvantage that, ultimately, only a subset of 
events (those that are near extremal) are used to generate the mass data, they have a number 
of obvious advantages. To begin with, the procedure is straightforward in its conception and 
implementation. Secondly, and perhaps most importantly, the procedure has little or no model- 
dependence. Were it to be established that a relatively pure sample of events from a particular decay 
channel could be isolated, then the subsequent analysis is pure kinematics. No assumption about 
spins, couplings or matrix elements is required. Finally, we remark that another supposed benefit 
of endpoint techniques is that they place less stringent demands on the degree to which the detector 
acceptance and backgrounds must be understood: the location of the step within a Heaviside-like 
function may be determined even if the height of the step function has been modulated by a 
smoothly varying unknown acceptance and added to a smoothly varying background of unknown 
shape. In a typical experiment, provided one ends up with a reasonable number of events that are 
extremal or nearly so, one should expect to be able to measure endpoints. 

We note that this application of endpoint techniques to mass determination at the LHC is 
not new. For example, the endpoint of the dilepton invariant mass distribution has long been 
suggested as a good observable for constraining the mass-squared differences in supersymmetric 
cascade decays [7]. Here, by contrast, we shall attempt to use the procedure to provide a direct 
measurement of the absolute mass scale, rather than differences in masses or differences in squared- 
masses. We do this by generalizing the method described above for the PF-boson. 

There are two obvious ways in which we would like to generalize the method used to determine 
the W mass. Both are motivated, in part, by our theoretical prejudice as to what we expect, or 
hope, to discover at the LHC, and both encounter an immediate obstruction. The types of events 
we will be considering are summarized in Fig. [TJ 

Firstly, we should like to generalize to cases where the invisible daughter has unknown, but 
non-negligible mass (Fig. f]Ta,c)). One theoretical prejudice for this is that any viable dark matter 
candidate ought to be electrically-neutral, colour-singlet and long-lived, hence invisible in the 
detector. Decays to this particle would be exactly of the type discussed above, except that the 
dark matter candidate, by definition, has non-negligible mass (perhaps of the order of the Fermi 
scale or greater, if the dark matter is both thermal and weakly-interacting). In such a case, even 
more kinematic information is lost, because both the momentum and the mass (and, ergo, the 
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FIG. 1: An illustration of the types of decay considered in this paper. We consider both single- and pair- 
production of a parent particle which can decay into one invisible daughter of unknown mass (represented 
by the lightly shaded arrows) and an arbitrary number of visible daughters. In the analysis and conclusions, 
we will find it necessary to draw a distinction between (a,b) events in which the parent (s) recoils against 
a system (large yellow arrow) with significant transverse momentum (SPT), and (c,d) events in which the 
parent(s) recoils against nothing and has zero transverse momentum (ZPT). The circular blobs serve to 
disguise the nature of the mechanism(s) by which each decay takes place, as our analysis is general in this 
regard. See Fig. [2] for some special cases that are examined in more detail in the text. 



energy) of the invisible particle are unknown. Hence, tut is no longer an observable. 

Secondly, we should like to generalize to cases where the parent particles are pair-produced 
(Fig. Pi7b,d)). This is the expectation in, for example, supersymmetric theories with conserved 
imparity, or in little Higgs models with T-parity. In such cases, each of the pair-produced parents 
of mass mo can produce an invisible daughter of unknown mass mi among its decay products. 
Should this happen, the transverse momenta of the invisible particles can no longer be individually 
inferred from the missing transverse momentum; only the sum of the transverse momenta is now 
constrained. It is, consequently, not possible to observe the transverse mass of either pair-produced 
particle. 

The generalization of tut to events with pair-production was first made in [8j |9], where an 



algorithm was given for constructing an observable, called mT2(x)i which is a function of the 
transverse momenta of the observed particles and the missing transverse momentum, and whose 
definition is based upon the individual (and unobservable) transverse masses. (An entirely different 
way of considering events of the same type was suggested in [TO].) Note that, in the construction of 
the rriT2{x) variable, it was necessary to introduce a parameter \ which represents a hypothesis for 
the true mass mi of the invisible particles generated in the decay of each primary parton. In that 
sense, m,T2{x) should strictly not be called an event variable, but rather an "event function". A 
given event generates not just one value of rxiT2- Instead it generates a distribution of tot2 values, 
one for each value \ taken as a hypothesis for the mass of the invisible final-state particles. When 
the hypothesis % is chosen to be equal to the true mass m\ of the invisible particles, the observable 
m T2{x I X = rn i) shares two properties enjoyed by rriT, namely (1) that it is bounded above by 
the mass mo of the decaying particle, and (2) that this bound may be reached in events in which 
certain constraints are met, such as the rapidity of each invisible particle matching the rapidity of 
the sum of its sister decay products. The upper endpoint over events (indicated by a circumflex) 
of m.T2, m T2(x I X = m l)> may thus be used as an estimator for tjiq. 

However, without knowledge of mi, what may be said about the dependence of m,T2{x) on 
the mass hypothesis %? Without knowledge of mi, is one forced to treat rh,T2{x) as providing a 
one-dimensional constraint in the (mo,mi)-plane parameterized by x- Is it the case that the m,T2 
variable constrains not mo and m\ but only a relationship between them? 

Cho et al. [lj considered "CASE 6v ZPT" and "CASE 6s ZPT", as defined in Fig. |J Their 
claim, based on simulations and examples, but made without proof, is that the rriT2(x) observable, 
considered as a function of the hypothetical mass of the invisible particle %, has a kink (that is, 
is continuous, but not differentiable) precisely at the point where the hypothetical mass equals 
the true mass. In [2 , one of us examined "CASE 2 SPT" analytically and proved the existence of 
a kink at x = rn i- I n Section III we analyze all the cases of Fig. [T] and Fig. [2] including, but 



not restricted to, those of [T] and [2]. We go on to produce example Monte Carlo distributions in 



Section IV including examples of CASE 6v ZPT and CASE 6s ZPT (Fig. 10 12) and CASE 2 SPT 



(Fig. l6b) which illustrate the sorts of kinks found in pQ and [2] respectively. 

III. ANALYSIS 

Our notation for a single particle decay is as shown in Fig. [2] The decaying parent is particle 
0, the invisible daughter is particle 1, and the system of visible daughters is labelled N. 




FIG. 2: This figure shows six particular forms of event which are used to illustrate points made throughout 
the analysis and in the conclusions. They are all special cases of the generic event-types shown in Fig. [I] 
Each case will come in both SPT and ZPT forms (i.e. with and without net transverse momentum as defined 
in Fig. nl although for clarity we omit the recoiling system from the diagram above. This figure also defines 
the numbering system used to denote particles in these special cases. The invisible final state particles are 
represented by a lighter colour, and are numbered either 1 or 6. 

We denote the energy of particle i by Ei, its transverse momentum by pi, the magnitude of 
its transverse momentum by pi, its longitudinal momentum by %, and its transverse energy by 
d- The mass of particle i is rrii and the transverse energy is defined in terms of the transverse 



momentum by e% = \ pj + mf. For the case of the invisible particle 1, we will sometimes need 



to use the hypothesized mass x an d the corresponding transverse energy, ei = \Jp\ + X 2 - The 
antisymmetrized product is denoted by -EVi^l = -£d<Z2 — Eiq\. Finally, we put a circumflex over a 
function which is evaluated at an extremum. 

A. Single Particle Production (CASE 2, CASE 3v, CASE 3s, etc.) 

Consider the decay of particle, 0, to an invisible particle, 1, and an iV-particle system of visible 
particles, % with individual mass m; and invariant mass mjv- The number of transverse dimensions 
is arbitrary. We first assert that tun S [to<, m>], where, for a true (N + l)-daughter point decay, 

m < = EiTrii, (7) 

m > = too — toi. (8) 



For decays of cascade type, which proceed through intermediate on-shell states, m < and m> will 
be larger and smaller, respectively. For a two-daughter decay, m< = m > = 77*2, whereas for a 
three-body cascade decay to massless visible daughters, occurring via an on-shell state of mass M, 

m< = 0, (9) 



/ (mg-M2)(M2-mj) 

m > - V ^2 • ( 10 ) 

For further details, see e.g. [71 ITT]. 

The problem is to extremize the function 

fix 2 ) = X 2 + m 2 N + 2{eie N - pi -p N ), (11) 

which is the square of the transverse mass one would calculate when assuming a hypothetical mass, 
X, for the invisible particle. 1 The extremization is subject to various constraints, namely conser- 
vation of energy-momentum, and the mass-shell conditions. We consider two distinct situations. 
In the former, the decaying particle can have arbitrary transverse momentum with respect to the 
laboratory frame (SPT). In the latter, the decaying particle has zero transverse momentum in the 
lab frame (ZPT). 

1. Single-particle SPT production fCASE 2 SPT, CASE 3v SPT, CASE 3s SPT, etc.) 

In the SPT case, a minimal set of constraints can be written as 

ml = m\ + m 2 N + 2(E 1 E N - pi • pat - <?i<?jv), (12) 

m\ = El-p{-<£, (13) 

m 2 N = E 2 N -p 2 N -q 2 N . (14) 



1 We could also choose to extremize the same function, but with the invariant mass mjv replaced by the transverse 
invariant mass of the visible system. But the transverse invariant mass is also valued in [m < ,m > ], and the extremal 
values of / defined in this way are the same. Moreover, there are two obvious disadvantages in using the transverse 
invariant mass. The first is that it only attains its boundary values, m<, when the relative rapidities of visible 
particles vanish. There are thus fewer events near the extrema for a given sample size. The second is that it does 
not use available information, namely the longitudinal momenta of the visible particles. 
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We introduce Lagrange multipliers Ao,i,jv and vary with respect to I£i,jv>Pi,iV>9i,iV and m^ to 
obtain the extremization equations 

= XqEn — XiEi, 

= XqEi — XnEn, 

= A <?jv - X x qi, 

= X qi — XxqN, 

= - P7V (l + Ao) + (Ai + ^)pi, 

ei 

= Pl (l + A ) + (X N + —)p N , 

e N 

= m A rfl + A + A A r + — V (15) 

\ eN t 

Note that we assume that none of e\, e± n vanish so that we can freely multiply or divide by any 
energy or transverse energy. Now, the first four of these equations together imply that 

X 2 = AiA^v, (16) 

and that either 

= A = Ai = X N , (17) 

or that 

= E [lQN] , (18) 

which, just as in Mb, implies that the relative rapidity of the visible system and the invisible particle 



should vanish. In the former case, equations (15) imply 



= eipjv - ejvPi, (19) 

or equivalently that 

mNPi = XPN- (20) 



But using (19) in (11) then yields 

f = (x + m N ) 2 . (21) 
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In this case also, the last of equations (15) has no solution for m^ G [ m <, my], which implies that 



there is no stationary point in rajy. Thus, the extrema of / must arise at the boundary values of 
rriN- The extremal values are 



fj = { x + m< ) 2 i 
fn = (x + m>) 2 - 



(22) 



For a true (./V+l)-daughter point decay, 



fi = X + ^im, 
fn = x-mi + m . 



(23) 



It is straightforward to check that all of equations (15), bar the last one, have been satisfied. One 



should also check that the constraints, (12), can be satisfied. They can, for some given value of p±. 



In the case m> = mo — mi, for example, one finds that 



2 2 

2 X m( 

Pi ~ 



(24) 



(x-™i) 2 ' 

The singular behaviour at x = m i arises, because there are extrema for all values of pi at this 
point. 



The second possibility ( 18 ) is that = E^g^i , in which case the first constraint in ( 12 ) simplifies 



to 



ml = m\ + m 2 N + 2(e 1 e N - pi • p^v)- 



Now the penultimate two equations of (15), together with (16), imply that 



= (eipAr - eArpi)(eArAo - eiAi). 



(25) 



(26) 



One of these two factors must therefore vanish; if it is the first, then we are led back to (19) and 
the extremum (21). If it is the second, we find, from the first equation of (15), 2 that 



E\eN = .EWei. 



(27) 



Now since = Enqm in this case, we can square and add {E\q^) 2 to the left-hand side and (E^qi) 



to the right-hand side to obtain 



ei = ei, 



(28) 



We assume all Lagrange multipliers non- vanishing, otherwise we end up with the extremum obtained previously. 
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or rather 



X = mi, 



(29) 



for which / = m^. So this possibility generates the well-known maxima of / at \ = 7ri ii with 
= E {1 q N] . 

We should also consider whether / takes extremal values elsewhere on the boundary. Indeed, 



we already saw in (22) the extrema occurring at the boundary values of the invariant mass mjq. 



There are two other boundaries that are consistent with the constraints, namely the boundary at 
large q\N and the boundary at large £>i,7v (both of which imply large Ein via the constraints). 
It is easy to see that the boundary at large q±^ does not give rise to new extrema for /, since / 
has no explicit dependence on qi^jy. Let us look instead at large values of the transverse momenta 
p\ and pn- We should still satisfy the extremization equations obtained by variation with respect 
to Ei^n and q\ n, and so we still find that either the Lagrange multipliers or the relative rapidity 



should vanish. In the former case, we reproduce the extremum (22). In the latter case, the first 



constraint in (12) becomes 



and / can be written as 



ml = m\ + m 2 N + 2(e 1 e N - pi • pjv). 



fix 2 ) = ml + x 2 ~ mj + 2eAr(ei - e{] 



At large p\ and pn, equation (30) becomes, at leading order, 



= PiPn ~ Pi • PN, 



(30) 



(31) 



(32) 



which is solved by choosing pi and pjy to be parallel. At next-to-leading order, (30) becomes 

2PN 



2 Pi 2 2 2 

mr- h m N — =m Q - m 1 - m N , 

Pi PN 



(33) 



whence 



PN 

Pi 



ml — ml — m 2 N ± 



Y (ml -ml- m%) 



2 ^2 _ 4 m 2 m 2^ 



2ml 



(34) 



To understand how the two values arise, note that in the rest frame of particle 0, the magnitudes 
of the momenta, pin, are completely fixed (for a given value of the invariant mass m^). The two 
values above are obtained by performing an infinite boost either parallel or anti-parallel to pi (in 
the rest frame of 0). 
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Now let us expand / in (31) to next-to-leading order. We find 



/-^HMA— <»f>(l+~) : 



m 2 Q + (x 2 -m\)(l + 



m\ — m\ — m 2 N ± * (ttiq — mf — m 2 N ) 2 — Am 2 



2^,2 



m 



N 



2mj 



(35) 



(36) 



Again, the extrema with respect to tun are obtained at the boundary, such that the extrema are 



rrin — mr — m; 



+ ^ 



i<) 2 — 4m 2 m 



2^2 
< 



2m? 



flll = m 2 + ( X 2 -m 2 )(l+'-^ 

f m 2, fv 2 ^ ^ i m o ~ m i ~ m > + vVo ~™(- ™j) 2 ~ 4mfm| 
7rv = m + (x -"T-iJUH — 2 



fv = m 2 + ( X 2 -m 2 ) l + ^> 



mn — mT — m: 



^/(toq — m 2 — m>) 2 — 4m 2 m 



2m? 



fvi = m 2 + ( X 2 -m 2 ) l + ^> 



m-n — mT — m; 



^/(toq — m? — m<) 2 — 4m 2 m 



2ml 



(37) 



There are thus four possible extremal values for / which occur at asymptotically large momenta. 3 
They obey the following order relations above the kink (the order relations are simply reversed 
below the kink) 



fin > fiv > fv > fvi- 



The central relation is trivial; we prove the other two by establishing that 



(38) 



ml 



m 1 — m N ± \/ { m o ~ m i ~ rn %) 2 ~ 4m\ 



2™2 



m 



N 



(39) 



is monotonically decreasing on m^ £ [0, oo] for the + branch, and monotonically increasing for the 



branch. Indeed, differentiating with respect to m 2 N yields 



-1=F 



rn + m l ~~ m Af 



J (ttiq — mf — m 2 N ) 2 — 4m 2 m 2 N 



(40) 



Furthermore, the magnitude of the fraction is greater than one for all values of m 2 N . Hence (39), 



has the claimed monotonicity properties, and the order relations in (38) follow. 



For a true point decay with massless visible daughters, the gradients are given by 1, mo /mi, 
mo/ mi and (mo /mi) 2 . 



In summary, the possible extrema for / are given by fi-vi in equations (22) and (37). 



Two in the case of a two-daughter decay with m < = m ; 
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Up until now, we have said nothing about the nature of the extrema of /, that is whether they 



correspond to maxima, minima, or saddle points. Now, it is clear from (11) that the extrema of 
the type 



fl,ll = (X + m^) 2 (41) 

correspond to minima in phase space. That is, if we shift the momenta slightly, subject to the 
constraints, then / will increase. Moreover, if we increase mjv from m<, then / will also increase. 
Thus 

// = (x + m<) 2 (42) 

is a minimum. On the other hand, if we decrease m^ from m>, then / will decrease, implying that 

A/ = (x + ™>) 2 (43) 



is a saddle point. What about the asymptotic extrema? Well, we know that the — branch in (35) 
corresponds to a minimum value for the coefficient of (x 2 — Wi), with respect to variations in the 
momenta, whereas the + branch corresponds to a maximum. But we also know that the coefficient 
is monotonically decreasing as a function of to^t on the + branch and monotonically increasing 
on the — branch. Putting this together, we see that, above the kink, fjjj is a maximum, fyj is 
a minimum, and fiy and fy are saddle points. In contrast, below the kink, fyj is the maximum 
and fin a minimum. 

Since there is only ever a single maximum, for any value of x, it must be the global maximum. 
Since there are always two minima, it is not immediately clear which is the global minimum. In 
the case of a point decay to massless visible particles, it is easy to show that // is always the global 
one. 

In summary, the global maximum for \ < m i 1S given by 

it i \ 2/2 2\ /-, m 2 , — m 2 — mi — Jim 2 , — m 2 — mi) 2 — Am 2 mi \ 
f(x\x<™ 1 ) = m 2 + ( X 2 -ml) (1 + -2 1 < v ° 2 x <> l -^ , (44) 

whereas for \ > m i> it is given by 



+ \/(ttIq — m 2 — m 2 ^) 2 — 4m 2 7 



2 _ ^,2 _ _2 , . fT—2 _^2 _ rrn 2\2 _ A^.2^2 



rrin — mt — mz- + \/ [mZ — m? — mz- r — 4m?m 



/(xlx>mi) = ^o + (x 2 -m?) (1 + -^ ^ V " 2 — ^^ • (45) 

These functions coincide at % = mi, but have different gradients there, giving rise to a kink. 
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2. Single-particle ZPT production fCASE 2 ZPT, CASE 3v ZPT, CASE 3s ZPT, etc.) 

To study decays in which the decaying particle has zero transverse momentum, we simply 
append the constraint 



= Pi + pjv, (46) 



to the set ( 12 ). As before we find that either all Lagrange multipliers vanish or the relative rapidity 



vanishes. Then, we find from (12) and (46) that 



/ 9 9 9\9 ^99 

(Mq — m\ — m N ) — Am\mj, 
Ami 



i>i - " ' Ef " " . (47) 



whence 



and 



Thus, 



m 2 + m\ - m 2 N 
61 = 2m7 (48) 



t2 _ 

2 mo 



eN= m -m 1 + m N _ ^ 



, 9/9 9 N ml-m\ + m 2 N ( „ / mX + mf - m 2 N \ 2 m3, + m\-m 2 N \ 

f = ml + {x 2 -m\) + ^ N \j^-m\) + { ^ " ) - ^ " j- 

(50) 
Note that the radicand is always positive-definite for x > 0, and < tojv < mo — m i (the largest 
possible endpoints for mjv), such that / is real. 

We have not yet considered the behaviour of / as we vary m^. Let us now prove that / is a 
monotonically increasing function of m^ for x > m ii an d a monotonically decreasing function for 
X < mi, such that the global maximum is obtained at the upper endpoint of m^ for x > m i an d 
at the lower endpoint for x < m i- To do so, it is convenient to define the quantities 



A = m 2 



o- 

B = 4m 2 ( X 2 -m 2 1 ), 
C = m\ — m 2 N , 
D = m 2 + ( X 2 -m 2 ), (51) 

which are all positive semi-definite apart from B > —Am^m 2 . In terms of these quantities, / may 
be written as 

A-C 



f = D + ^ r ^B + (A + C) 2 -(A + C)), (52) 
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such that 

df 



B + 2C(A + C) - 2CyjB + {A + C) 2 ) . (53) 



dm 2 N 4A^B + {A + C) 2 
To exhibit the claimed monotonicity properties of /, we wish to show that, for B ^ 0, 

B + 2C{A + C) - 2C^B + {A + C) 2 ) § 0. 



Now, since AC and (A + C) 2 are both greater than Am^m 2 , we have that B + 4AC > and that 
B + (A + C) 2 > 0. From the first of these, B «g implies 

£( J B + 4AC) >0, (54) 

which itself implies, after adding an identical term to both sides, that 

{B + 2C{A + C)) 2 ^4C 2 {B + {A + C) 2 ). (55) 

Since C > and B + (A + C) 2 > 0, we can safely take the square root on the right-hand side to 
obtain 

\B + 2C{A + C)\§ 2C^B + {A + C) 2 . (56) 

Now B + 2C(A + C) can be either positive or negative. If the latter, then B is necessarily negative, 
and, moreover, B + 2C(A + C) — 2C\J B + (A + C) 2 is the sum of two negative-definite terms and 
is negative. If the former, then (|56|) becomes 



5 + 2C(4 + C) - 2CyjB + {A + C) 2 ^ 0. (57) 



In either case, we have the result claimed, namely that B ^ implies I B + 2C(A + C) — 
2Cy ' B + (j4 + C) 2 I 5s 0, which in turn implies that / is monotonic, either increasing or decreasing 
as specified above. Thus, the global maximum of / for \ < m\ is given by 

f(x\x < m i) = 



m 2 , — m 2 + m 2 ^ ( I „ 2 / to 2 , + m 2 — m\ \ 2 itiq -\- m 2 



m 



2 



whereas for x > mi, the global maximum is given by 
/(xlx > mi) = 



^o + (x 2 -m 2 ) + * J (x 2_ m 2 ) + ( i >; ) - °i 



2mo \ V ^ 2mo ' 2mo 

(59) 
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These two values coincide at x = rn i but have different gradients there, iff. m > ^ m < . This 
condition is trivially broken in the case of decay to a single visible daughter, in which case m> = 
m < = m<i. In all other cases, we will get a kink. 

It is worth observing that the maxima of the ZPT case never coincide with maxima of the 
SPT case, not even as local (rather than global) maxima. Indeed, the only maximum of the ZPT 
case which coincides with an extremum of the SPT case occurs for a true point decay, for which 



rriN = mo — mi and pi = 0. But, as we pointed out in the previous subsubsection (IIIA1), this 
extremum is a saddle point in the SPT case. 

B. Further analysis of single particle decays (CASE 2, CASE 3v, CASE 3s, etc.) 



Section III A| provides a comprehensive description of the extremal values of / for single particle 



decays, but is not terribly intuitive. 

In order to gain some intuition as to the physical origin of the kinks and extrema, in this 
subsection we investigate the dynamics by first looking at the decay kinematics in the rest frame 
of the parent particle, applying the physical constraints there, and then boosting the system of 
daughters to the lab frame. As before, the parent particle, '0', decays to an invisible particle, '1', 
and a system of visible particles, i N\ (En, p^, gjv) can be either the four-momentum of a single 
particle in the two-daughter decay case, or the sum of four-momenta of all the visible daughters in 
the iV-l-l-daughter decay case, with N + 1 > 2. 

We have shown already that maxima of / occur when the visible and invisible systems are 



produced with vanishing relative rapidity (18). Since it is the maxima that we are interested in, 
for the rest of this section we restrict ourselves to vanishing relative rapidity, in which case the 
transverse mass-squared is given by 

f = m 2 + (x 2 -m 2 l )+2e N (e 1 -e 1 ). (60) 

Moreover, since we are interested in the existence of kinks near x = m\, we simplify subsequent 
expressions by expanding / as a Taylor series in x 2 i about x 2 = m i- Thus 

f = m l + {x 2 -ml)(l+ e -£}+... (61) 

We remind ourselves that since we are interested in the maximal value of /, for values of x 
greater than mi, we must take the maximal value of 

«^1 + — , (62) 

ei 
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whereas for \ l ess than mi, we must take the minimal value of a. There will be a kink in / if the 
maximal and minimal values of a are different. 

We first consider the decay in the rest frame of the parent. In that frame (denoted by primed 
quantities), for fixed values of mo, m\ and tun, the momenta and transverse energies of the daugh- 
ters are fixed: 

{P) = |Pl > 2 ' " 4m§ ' (63) 

&? = &? + ml (64) 



We now boost back to the lab frame. We define (3 = |po/ e o| to be the transverse velocity of 
the parent particle in the lab frame. We define 6 to be the angle, as measured in the rest frame 
of the parent, between the velocity vector of the lab frame and the transverse momentum vector 
of the visible particle system, p^r. This means that, to an observer in the lab frame, the condition 
6 = corresponds to the invisible particle being thrown forwards in the decay, i. e. in the same 
direction as the parent's velocity, while 9 = ir corresponds to the invisible particle being thrown 
backwards, i. e. against the parent's velocity. 

A simple Lorentz transformation gives the lab-frame transverse energies and momenta of the 
daughters, 

eN = l(e' N -p'/3 cos 9), (65) 

ei = r y(e 1 + p j3 cos 9) , (66) 



i 



where 7 = (1 — (3 2 )~^ is the Lorentz factor associated with the velocity (3. This means that a is 

given by 

e' H — p'B cos 9 

a = 1 + — — — 

e i + P'fi cos 9 

2ml 



m o + m i ~ m N + 2moP'/3 cos 9 ' 

(67) 



2ml 



ml + m\ — m 2 N + [3 cos 9 J (m 2 , — (mi — m^v) 2 ) (m 2 , — {m\ + mjv) 2 ) 
Equation 67 tells us everything about the gradient around x 2 = m i f° r an y two-daughter, 



three-daughter, or indeed iV+l-daughter decay. We recall that N represents the system of visible 
daughters and has a mass equal to the invariant mass of the visible daughters. For the iV+1- 
daughter decay, with N + 1 > 2, the visible-daughter invariant mass, tojv, can take any value 
within the range rriN £ [m < ,m > ]. We recall that, for a true point decay, m< is equal to the sum 
of the masses of the visible daughters produced in the decay, and m> is mo — mi. 
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We can see that, if some finite experimental sample contains events in which either the parent 
is highly boosted, such that (5cos6 has values significantly different from zero, or itin spans some 
reasonable range, then the maximal and minimal values of the gradient obtained from the sample 
will be significantly different, and an observable kink in / can be expected at x = rn i- Indeed, in 
principle, we need just two events with vanishing relative rapidity to generate a kink at x = m i- 

Thus, the presence or absence of the kink at (mo, mi) is not, in fact, contingent upon whether 
our experimental sample contains events that correspond to the global maxima, whether SPT or 
ZPT. Rather, we simply need a reasonable number of events with vanishing (or near-vanishing) 
relative rapidity and significantly differing values of either m^ or (5 cos 6. Either of these will suffice 
to generate a kink, and indeed our later simulations will illustrate this rather well. 

That said, it is important to remark that, if our sample does not contain events at or close to 
the global maximum, then it is possible that further kinks will arise at x ¥" m i- The reasoning for 
this is as follows. Suppose our sample contains events with vanishing relative rapidity, giving rise 
to a kink at x = m i- If an event corresponding to the global maximum is not contained in the 
sample, then it is possible that other events in the sample, with non-vanishing relative rapidity, 
can exceed the maximum value of / generated by the events with vanishing relative rapidity. Now, 
this cannot occur at % = mi, because we have saturated the global maximum at this point. So, 
the kink at (mo, mi) cannot be erased in this way. But the events with non- vanishing rapidity can 
exceed the maximum value of / generated by the events with vanishing relative rapidity at values 
of x 7^ m i- If this does indeed occur, then a second kink will be generated at the value of x ¥" m i 
where the values of / coming from the two different types of event coincide. 

So it is possible, albeit unlikely, that a finite experimental sample of events will give rise to 
spurious kinks, leading to a discrete ambiguity in the extracted values of the masses mo and mi. 
We note, however, that this phenomenon is not observed in any of our simulations. 

Let us put this issue to one side, and return to our discussion of the kink at (mo, mi). In the case 

where the range of values of rriN is small — and in particular CASE 2 where m^ is single-valued 

- we must rely on different values of (5 cos 9 to produce a kink. The global-maximum gradient on 

either side of the kink is obtained asymptotically as the transverse momenta become both large in 

magnitude and collinear. 

Some other important properties of the extremal gradients are as follows. 

• The maximum and minimum gradients will occur when (5 cos 6 — * =pl. The values of these 
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/3cos 




400 



FIG. 3: The surface shows the gradient, a = j^, near \ 2 = m i as a function of /3 cos 9, and of m/y, the 
invariant mass of the visible particles produced in the decay. The values of the other masses are mo = 500 
and rti\ — 100. 



extremal gradients are (c/. Eqs. 44 and 45) 



o.~ 



1 + 



m 



±V^ 



1 - u /( rrin — m 2 



m 



2 < ) 2 - Am\ 



2 m1 



2m\ 



In the limit where all the visible particles are massless, m< = 0, these become: 



a 



a 



raoY 
mi/ 



1. 



(68) 

(69) 
(70) 



• When m> = itiq — mi, then p' = 0, and the invisible daughter is at rest in the parent's 
rest frame. The gradient is then independent of the motion of the parent (parameterized by 
(5 and 6) and takes value mo/mi. In CASE 2 this configuration is only obtained if masses 
are precisely at the threshold. In the N + 1-daughter point-decay case with N > 2 (e.g. 
CASE 3v) , it represents the kinematic configuration in which the visible daughters have their 
maximum possible invariant mass. This configuration cannot be realised in cascade decays, 
such as CASE 3s. 

In Fig. [3] we show the gradient as a function of (3 cos 9 and mjy, for the particular values of 
mo = 500 and mi = 100. As expected, the maximal value of the gradient occurs at (3 cos 9 — ► — 1, 
rriM = 0, and is rrig/mf = 25. The minimal value occurs at /3cos 9 — > +1, mjy = 0, and is 1. When 
rriTv = fno — m i = 400 we obtain gradient mo/mi = 5 for any value of (3cos9. 

Note that in a real collider experiment there will usually be very few events near (3cosQ = ±1 
since these events require /3 — ► 1 and hence asymptotically high energy. Indeed, for the common case 
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of heavy particles being produced near threshold, most events will have (3 ~ 0. These threshold- 
production (ZPT) events have gradients 

<* p=0 = l + ^, ( 71 ) 



2m\ 



m + m l ~~ m 7V 



(72) 



For the iV+1-daughter point decay with m< = Ejmj and m > = tuq — mi, the ZPT extremal 
gradients are 

„max "*0 ( , 'tq\ 

^=0 = — , (73) 

°n = m 2 + Tf - m 2 < • (74) 

Returning to the SPT case (in which the global maximum is obtained at asymptotically large 
momenta) we should like to understand how close to the global maximum we can get in events 
with large, but finite, momenta, such as we might hope to achieve in an experiment. To this end, 
let us ask what the relative decrease is in the gradient of / that results from a small shift in /3 cos 9 
away from the maximal values of ±1. For (5 cos 6 = 1 — <5+, we find that the relative change in the 
gradient is given by 

whereas for (3 cos 9 = — 1 + <L. , the relative change in the gradient is given by 

4^-,- (76) 

e\ — p 

Thus, the relative change for a given shift in the 5s is always greater above the kink than it is 
below, irrespective of the mass values. So in a finite sample of events, in which the distribution of 
cos 9 is roughly uniform, we always expect to get closer to the maximum below the kink than we 
do above it. 

A graphical representation of this result for particular mass values can be seen in Fig. |4j where 
rriT = \J f is plotted as a function of cos0 and (3 for two particular values of %. Firstly, we see 
that the maxima are obtained at the expected values of 9 and [3. We also see that, for angles close 
to the respective maxima, the fall-off is much greater above the kink (x > mi, Fig. Ep) than it 
is below (x < m ii Fig. Ek). Indeed, a large fraction of the surface lies near the upper limit for 
X < mi, while only a very small fraction of events, with (3cos9 ~ —1, lie near the maximum value 
for x > m i- 
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(a) 



(b) 



FIG. 4: The surface shows the maximum transverse mass that would be found for a two-daughter decay 
when the hypothesized mass, %, of the invisible particle is (a) less than and (b) greater than its true mass. 
The invariant mass is plotted as a function of the transverse velocity, j3 — po/eo, of the parent particle 
in the lab frame, and of the cosine of the angle, 8, as described in the text. Decays in which the visible 
and invisible daughters are produced with vanishing relative rapidity will lie on the surface, while all other 
events will lie below it. In both cases the masses are: parent ttiq — 500; visible daughter 771,2 = 0; invisible 
daughter mi = 100. The hypothesized mass of the invisible daughter (x) is 50 in (a) and 150 in (b). 



We might therefore expect that in a finite sample of events, we are likely measure a maximal 
gradient close to the predicted one when x < m i- However we are likely to underestimate the 
maximum gradient when x > rn i because of the scarcity of events in this kinematic region. Our 



Monte Carlo simulations (described in Section IV) confirm these expectations. 



C. Pair production at SPT and ZPT (CASE 4, CASE 6v, CASE 6s, etc.) 



Now let us consider the case, relevant to supersymmetric theories with conserved R-parity, 
in which parent particles are pair-produced, and each parent decays into an invisible daughter 
particle and one or more visible daughter particles. This case is complicated by the fact that 
there are now two invisible particles, and thus two transverse momenta go unmeasured. Only 
their sum can be inferred from measurement of the total missing transverse momentum. Thus, the 
individual transverse masses for each decay are no longer observables, and we need to construct a 
new observable. 

A suitable observable was put forward in [HI [9], and is defined as follows. We label the pair of 
parent particles as and 5 and their respective invisible daughters by 1 and 6 (see figure p| . There 
are N and M visible daughters, respectively, with invariant masses tjin and thm- For simplicity, 
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we assume that the produced pair of particles have identical mass mo, that the invisible particles 
have identical mass mi. 4 

Following OS], we construct two transverse mass functions, namely 

f = x 2 + m 2 N + 2(e 1 e N -p 1 -p N ), (77) 

g = x 2 + m 2 M + 2(e 6 e M - Pe • Pa/ )• (78) 

Even if the mass % were known, these would not be observables, since one cannot separately 
determine the transverse momenta pi and p6 of the invisible particles in a collider experiment. 
One can only determine their sum, which equals the total missing transverse momentum, p, say. 
An observable can then be constructed as follows. Consider all possible partitions of the measured 
missing transverse momentum p into the unmeasured invisible transverse momenta. Only one of 
these partitions will, of course, be the correct one. For each partition, select the largest value of / 
and g, max(f,g). Finally, minimize max(/, g) over all possible partitions. 

In this way, one obtains an observable function of the hypothesized invisible mass x- This 
observable, called my 2 (x) has the same property as the usual transverse mass observable: at 
X = mi, it is bounded above by m 2 , [U [9]. 

To study m^ 2 away from the point x = m i> let us define the hypothesized invisible momenta 
of particles 1 and 6 in the trial partition to be p^ and pg. Since p^ + pg = p we can write 



/ = X 2 + m\ + 2J{p\f + X 2 Jp 2 N +m%- 2p* • p N , (79) 



9 = X + m 2 M + 2y (p - Pi) 2 + X 2 V Pm + m M ~ 2 (P " Pi) • PA/- (80) 

The prescription for constructing m 2 -^ given above instructs us to find the minimum with respect 
to variations of p*, of the maximum of /, g. There are three ways in which this minimum can arise 
(see Fig. ^J . Either (a) it is a minimum of one of the /, g that lies above the other one of g, /, or 
(b) it is a point at which / = g, or (c) it occurs at a boundary. In the case at hand, it is not difficult 
to show that /, g, considered as functions of p^, are unbounded above as p^ — ► ±oo and each has 
a unique minimum with /, g = (% + m,jv,A/) 2 - For /, this minimum occurs at XPN = ^-TVPi, and 
for g, it occurs at XPM = iti>m (p — Pi)- F° r minima of type (a) the value of rriT2 that results is 
given by x + ^iV,M, and thus the largest value of rriT2 that can arise is given by the largest upper 
endpoint of the two rriN,M intervals. Such extrema always lie below mo at x = rn i (except in the 
special case of a point decay, in which the upper endpoint, tuq — mi, gives rriT2 = mo). Thus, they 
do not, in general affect the behaviour at (mo, mi) in which we are principally interested. 



This is, for example, the case for supersymmetric decays to the lightest stable superpartner. 
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(a) 


X 

(b) 




FIG. 5: A diagram demonstrating that the minimisation over some parameter of the maximum of two 
diffcrcntiablc functions may occur cither at (a) a minimum value of one of them, or (b) when they are 
equal, or (c) at the boundary of the domain. 



We are thus left to consider minima of the type (b), for which / = g. Now, to find the partition 
of p that minimizes / = g for a given event requires us to solve f = g, which is a quartic equation 
in one of the assigned invisible momenta. Rather than do this, we note that the entire procedure 
of computing mj<2 by minimizing over partitions and then maximizing rriT2 over all events, can 
be phrased as a single extremization problem. Indeed, we wish to find the extremum which is 
a minimum of / (or g, it does not matter which) with respect to variations of the hypothesized 
transverse momenta and a maximum with respect to the true event variables, subject to the usual 
energy-momentum constraints, supplemented with the constraint / = g. 

Explicitly, we wish to extremize 



/ = rn% + x 2 + ^{e\e N - p} • Pat), 
where e\ = y/p* 2 + X 2 , subject to 

= m\ + x 2 + 2(e£ejv - Pi ■ Ptv) - m 2 M + X 2 + 1(e%e M - Pe ■ Pm) 
= -ml + m\ + m 2 N + 2{EiE N - pi • Pat - 9ig.iv), 

= -ml + m\ + m 2 M + 2(E 6 E M - P6 • Pm - <?6<?m), 

= -P + Pl +P6, 
= "P + Pi+P6- 

In the case of ZPT decays, we must further supplement these with the constraint 

= p + PAT +PM- 



(81) 



(82) 



(83) 



As always, variation with respect to E and q will force either the relevant Lagrange multipliers to 
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vanish, or the relative rapidities to vanish. That is 

= E^q N ], 
= E [6lM] ■ 



(84) 



We are interested in the latter case, since this will give rise to values of rriT2{x) which go through 
the point (mo, mi). In this case, the two energy-momentum constraints simplify to 



= -ml + m\ + m N + 2(eiejv - Pi ■ Pat), 
= -ml + ml + m 2 M + 2(e G e M - P6 ■ Pm)- 



(85) 



Let us now explicitly introduce Lagrange multipliers Mi, 2, 3 for the first three constraints in (82). 



The final two constraints we use to eliminate P6 and pg. We also add a vector Lagrange multiplier 



/i4 for constraint equation (83) in the case of ZPT. The extremization equations, obtained by 



varying with respect to p\, pi, Pat, Pm,P, t^n and tum are 



Di 

1 + ^i)[-^e N -p N 



'Pi — P 

Mil — ^ — e A f + Pm ), 



= jU 2 — e N - pjv J + A«3 e M + Pm 

V e\ I V cq 

,. , ,(PN~* *\ . (PN \ . 

1 + m) e 1 — Pi ) + M2 ( ei - pi J + £i 4 , 

-//1 1 e 6 - (p - Pi) I + M3 ( e 6 - (p - pi) + /x 4 , 

\e M > \e M 

-Mi — =i — ejv/ - pm + M3 eM -Pm + M4, 








eg " / ' "V e 6 

= (1 + m)[e\ + e N ) +M2(ei + e^v 
= -Mi (eg + e M ) + M3 ( e 6 + e M J , 

where it is to be understood that e§ and eg depend implicitly upon pi, pf and p. 



(86) 



We have not attempted to find the general set of solutions to (86) and (82). However, it is 



possible to show that, in the ZPT case, a solution of (82), (83), and all but the last two equations 
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of (86) is given by 



Pi = Pi, 

Pe =Pi> 

P6 = Pi, 

Ptv = -Pi, 



Pm = -Pi, 
1 
~2' 



Ml 

M2 =^3 
/i 4 = pi 



e l e 7\r e l 



ei2(eAr + ei) 2 ' 
(gi + ejy)(ei +ei) 
2ei(ei + eAr) 



7T17V = "T-Af • 



(87) 



In the above, pi is then fixed by the second of equations (86). Now, these values do not represent 



a stationary point of mx2, because they do not solve the last two equations of (86). If however, the 



two decay channels are identical, such that the endpoints of the m^ M intervals coincide, then we 
see that we obtain extrema at the upper and lower endpoints of the common invariant mass interval 
[m < ,m > ]. These extrema correspond to events in which both of the decaying parents are at rest 
in the lab frame, and in which each of the visible daughter systems obtains either its maximum or 
its minimum value. They are, thus, simply pairwise copies of the extrema we found for the single 



particle ZPT decays in subsubsection III A 2 They are, moreover, the extrema exhibited by Cho 



et al. in [T]. Since they are just pairwise copies of the extrema we found in subsubsection III A 2 
it is trivial to see that the resulting values of rriT2 generate a kink at (mo, mi) with m^ 2 given by 



m 2 T2 (x < mi) 



ml + (x 2 



and 



m\) + 



m- r 



m\ + m?^ 



2mo 




m% 



+ 



TTln 



mt 



mt 



2 mo 



ttiq + m 



m; 



2m,o 



ml 2 (x > mi) 



ml + (x 2 



rnf + 



m 



m\ + m> 



2m-o 



(X 2 ~ mf) + 



rriQ + m 



m: 



2 mo 



m.Q + m 



m: 



2m,Q 



(89) 
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Note that we have not shown that the extrema that we found correspond to the global maximum. 



To do so would require a complete solution of equations (86) and (82). But we can show that a 
kink is necessarily present nevertheless. Indeed, suppose there exists a maximum whose value of 
m,T2 is greater or equal to the values found above for some values of \. At x = m ii the rriT2 value 
of such a new maximum is necessarily mo, since we know that mx2(x = m i) 1S bounded above by 
mo- Now, considering values of x m the vicinity of x = m i> we see that any new maximum can 
only accentuate the kink at (mo, mi). This proves that rriT2 has a kink at (mo, mi) in the ZPT 
case, and indeed in any class of events that includes the ZPT such as a subset. Since the SPT case 
is such a class, we have also proven, as a corollary, that there is necessarily a kink in the SPT case 
as well. 

Note that the sole exception to this result is CASE 4, for either ZPT or SPT. In this case 
m < = m> = m,2 and the events described above do not generate a kink. This is just as was the 
case for the single particle ZPT case, CASE 2. Of course, this does not prove that there is no kink 
in this case, since we have not excluded the possibility of a higher maximum. Nevertheless, our 
later simulations will suggest that there is no kink in the ZPT case. 

Another immediate result that follows is that, just as for single parent decays, the extrema 
of the ZPT case we have identified are not extrema of the SPT case. Indeed, in the SPT case, 



we should solve equations (86) and (82), but without (83) and with /X4 = 0. But \i± = is not 



consistent with the solution we found in (87). Thus, configurations of this type are not extrema of 
the SPT case. 

What then, is the maximum of mn in the SPT case? To answer this, it is first useful to observe 
that the value of m^ 2 measured in an event cannot exceed the maximum of the individual transverse 
masses squared /, g that would be measured in the same event, if the invisible momenta could be 
determined. This follows from the definition of mx2'- the value of m^ 2 that one computes for a given 
event, viz. the minimum with respect to partitions of the invisible momenta of max(f,g), cannot 
exceed the value of max(/, g) that would be obtained with the true values of the invisible momenta, 
because one of the partitions must correspond to the true values of the invisible momenta. So m^ 
cannot exceed max(/, g) in a given event. This, of course, was the original motivation for the 
definition of rriT2 at x = mi, but it is true for all values of x- 

It follows that the maximal value of m^ 2 , obtained by considering some set of events, is bounded 
above by the maximal value of max(/, g) for the same set of events. If, furthermore, we can show 
that the bound is saturated, then we can infer the maximal value of m^ 2 for the set of events from 
the maximal values of / or g determined in the previous section. 
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This argument allows us to determine the maximal value of rriT2 i n the SPT case, where the 
initial transverse momentum of the pair is allowed to be arbitrary. Indeed, consider the event 
configuration in which each of the two decays in the pair corresponds to the global maximum of 



the single particle SPT decay found in Section III A 1 with arbitrarily large transverse momentum. 



Thus, the pair of decaying particles have equal and arbitrarily large transverse momentum, as 
measured in the lab frame. Furthermore, they both decay in the same way in their rest frames: 
for x < m i the invisible daughters are emitted in the forwards direction, and the visible daughter 
systems of invariant mass tti^t = rnu = r n< are emitted in the backwards direction. For x > thi, 
the invisible particles are emitted in the backwards direction in the rest frame. 

It is straightforward to show that the process of minimizing over the partitions of the invisi- 
ble momenta in this configuration singles out the partition with the true values for the invisible 
momenta. This shows that this configuration saturates the bound on vtiTi- 

In summary then, we have shown that in the case of identical pair decays, either SPT or ZPT, 
there is always a kink, except perhaps in CASE 4 ZPT. In the case of SPT, the extrema correspond 
to pairwise copies of the global maxima in the single-decay cases, and are the global maximum. In 
the case of ZPT, the extrema correspond to pairwise copies of the global maxima in the single-decay 
cases, but it is not clear whether or not they are the global maximum. 5 

IV. MONTE CARLO SIMULATIONS 

As a concrete illustration of this method, we perform simulations of events for particular sce- 
narios, analysing both single-parent cases, and also examples of pair-production of parent particles, 
as shown in Fig. [2j 

We generate Monte Carlo events by two different methods: 

• Firstly, we use the HERWIG[12, l2J|]3| Monte Carlo generator, with LHC beam conditions, to 
produce unweighted supersymmetric particle pair-production events. These events have the 
advantage that the transverse boost of the particles should bear some resemblance to that 
expected for real events, and so we use this generator when investigating SPT events with 
a "physically reasonable" px distribution. 6 We always generate the SUSY particles in pairs, 



-J 



It is worth remarking that the ZPT configuration in which the parents have arbitrarily large momenta, but are 
produced back-to-back, does not give a higher maximum. Indeed, the minimization over partitions in this case 
picks OUt 1TIT1 — y/f — X + m g • 
6 We recognise that it will somewhat underestimate the high-pT tail, as extra radiation ought to have matrix element 
corrections which will not be well modelled by HERWIG's parton shower. 
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but when we are considering the kinematics of single particle decays, we treat each decay 
independently, ignoring any missing transverse momentum with its origin in the other 'side' 
of the event. 

Secondly, to investigate the behaviour of the transverse kinematic endpoints under some 
unphysical but nevertheless highly instructive limiting cases, we use a home-grown toy Monte 
Carlo (MC) event generator. This toy MC allows us to either (a) generate events with the 
primary particle(s) having exactly zero net transverse momentum, thereby allowing us to 
examine the ZPT limiting case, or (b) generate events where the transverse momenta are 
non-zero (SPT) and entirely under our control. In the latter case, except where stated 
otherwise, we choose the pr distribution to have an extremely long unphysical tail (denoted 
"SPT large pt" in figures) in order to highlight the effects of large transverse momenta on the 
kinematic endpoint distributions. The toy MC generates all decays according to phase space 
alone, implements no showering and keeps all particles on mass shell. It defines an arbitrary 
mass scale 7775 = mo + m\ in terms of which its generation scheme is fixed as follows. In the 
single-particle production cases with ZPT, the parent particle starts its life at rest in the 
lab frame. In the ZPT pair-production cases, particles and 5 are generated back to back 
from the isotropic decay of an initial state at rest in the laboratory frame with centre-of- 
mass energy tuq + 777,5 + T, where the kinetic energy T differs from event-to-event, being an 
exponentially distributed random variable with mean 7775. In the single-particle production 
cases with SPT, the parent particle and an associated state with mass 7tt,jsr, against which 
it recoils, are produced back-to-back from an initial state at rest in the lab-frame which has 
energy 777,0 + T. Both ?77isr and the kinetic energy T differ from event to event, being random 
variables: ttt-isr is distributed exponentially with mean 777,5, while the kinetic energy T is 
distributed as the absolute value of a Cauchy-distributed 7 random variable with mean and 
scale 777,5. For pair-production with SPT, the two parent particles and an associated state 
with mass tojsr against which both jointly recoil, are produced in the isotropic three-body 
decay of an initial state at rest in the lab frame which has energy 777,0 + m b + T, where ttt-isr 
and T have the same distributions as they had for single-particle SPT production. 



7 There are two isolated cases in the paper where we replace this Cauchy distribution with an Exponential distribution 
with mean ms, thus significantly reducing the high pr tail to a more physically reasonable size. We do this on 
the two occasions in the paper in which we wish to use the toy MC to produce a less extreme and more physically 
relevant SPT distribution which is intermediate between its usual ZPT and usual SPT output. These two isolated 
cases are labelled "SPT (medium pT)" where they occur. 
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A. Single production with two-daughter decay (CASE 2) 
The predicted gradients at the kink are given for two-daughter decays as a function of the 



kinematic variables of the decay in (68 ). We recall that we expect qualitatively different behaviour 
around % = m i depending on whether or not the parent has transverse momentum in the lab 
frame. To be concrete: we expect to see a kink in the m^ ax (x) curve when the parents can have 
transverse momentum but do not expect to see a kink when parents are only produced at rest. 

We generated HERWIG events with slepton (£) pair production, followed by the decay £ — > £xi- 
The slepton mass was 300 GeV, and the Xi ma ss was 100 GeV. The slepton width was set to be 
much less than 1 GeV. The visible lepton is sufficiently light that it can be assumed massless in 
this context. We show in Fig. pk several distributions of rnr(x) = Vfix) f° r different values of x- 
We note that the upper kinematic limit is well saturated, making an experimental determination 
of the endpoint a reasonable expectation for this case. 

In Fig. [HJd we show the same information in a different format, now plotting the number of 
events (colour scale) with particular values of mr(x) (y-axis) as a function of x (x-axis). Using 
these distributions, we estimate the value of the maximum value, y/f(x) from the mid-point of the 
uppermost non-zero bin, which we denote m™ ax (x)- A graph of m™ ax (x) as a function of x shows 
clearly the kink at at (x = mi, my = mo); this is shown as a dashed line superimposed on Fig. rob. 



The gradients — — - 



- — were estimated from short line fits to that graph below and above 

the kink, and yielded gradients of 0.325 below and 2.56 above the kink. The predicted extremal 



values calculated for the example masses, using (|68|) and converted to 

dVf 



dx 



X df 



m dx 2 

X=m ±e 



(90) 

X=m ±e 



are ~ for x < rn i and 3 for x > m i- As discussed in Section IIIB we expect the measured value 
to be closer to the predicted one below the kink as compared to above it, since the high energies 
required mean that a much smaller proportion of events lie near the global maximum above the 
kink, as compared to below it. 

Note that the value of m™ ax we have used for the gradient determination is generated from 
the single maximal event at each value of x- I n an experiment which observes significantly fewer 
similarly-distributed events than we have generated (10 4 in Fig. pa-b), one would be unlikely to 
find an event at this limit. A typical experimenter observing, say, a factor of ten fewer events 
without background (or for whom e.g. ten events are required to get a statistically significant 
determination of the position of the end-point due to the presence of backgrounds) might expect 
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FIG. 6: (CASE 2) (a) Example distributions of the transverse mass, mxix) — Vf(x): f° r the two-daughter 
decay with mo — 300, mi = 100, and m-i = (HERWIG). The different histograms represent the same events 
but have different values of x- The thick solid line is for x — m ii while the thinner broken lines arc for other 
values of x uniformly spanning the range < x < 200. (b-d) The plots show the distribution of the number 
of events (colour scale) which produce particular values of m*r(x) — \/f(x) (y _ax is) for different values of x 
(a;-axis). (b) HERWIG events. The upper limit, m^ ax (x) is indicated as a dashed line. The short solid lines 
are fits to m™ ax (x) above and below x = m i- ( c ) F° r the phase space Monte Carlo in which the mother 
has been constrained to be at rest in the lab frame (ZPT). (d) For the phase space Monte Carlo in which 
the mother can have large transverse momentum (SPT) in the lab frame. 



to observe a maximum a little below our dotted within the neighbouring cyan regions of the plot. 

If we restrict the parent to have no transverse momentum in the lab frame (Fig. pc, ZPT) we 

find that the kink disappears. In its place we see a smooth curve passing through the point where 

the kink used to be, with the gradient at that point measured to be ~ 0.60. A plot where the 
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parent's transverse momentum, pr, can take very large values (Fig. pp, SPT) shows a slightly 
more prominent kink than the HERWIG plot. As we increase the parent's pr (ZPT— >HERWIG^SPT) 
the fraction of events contributing to the kink in a significant way increases. Both the absence of 
the kink for CASE 2 at ZPT and the increasing number of events found near the kinematic limit 



(Vf) as Pt (and hence (3) increases confirm the predictions of (67). The gradients measured on 
either side of the kink in Fig. [6]1 are found to be 0.335 and 2.99. These gradients are closer to the 
asymptotically achievable values of 1/3 and 3 than the gradients seen in the HERWIG events, in line 
with expectation. 

B. Single production with point-like three-daughter decay (CASE 3v) 
For the three-daughter decay, the range of allowed values of m^ means that we expect to see a 



kink regardless of whether the parent has any transverse momentum in the lab frame (67). 



As for Section IV A| we generate events with HERWIG, but this time for the initial process of 



gluino pair production. The gluinos then decay according to g — ► qq~Xi- The mass of the gluino 
was 500 GeV, and that of the Xi> 100 GeV. Again, the quarks can be considered massless for our 
purposes, since for simplicity we only allow decays to light quarks. 

In Fig. V7r we plot distributions of mj- = y/f for this three-daughter decay for different values of 
X- The upper limit is well saturated, and the edge of this distribution, m™ ax (x) also shows clearly 
the kink at {x = rni, my = mo) (Fig. 7p). The fitted gradients — ^ — are 0.223 below and 2.1 
above the kink. The predicted gradients are 1/5 below the kink and 5 above the kink for SPT 
events, and 5/13 ~ 0.3846 and 1, respectively, for ZPT events. The observed gradient on either 
side of the kink lies between the SPT global maximum and the ZPT maximum. This reflects the 
fact that the gluino has finite, but non-zero, pt in the lab frame. As expected, we get closer to the 
SPT maximum below the kink than we do above it. 

As predicted by the analysis, for the three-daughter case, the kink is retained in ZPT events 
where the parent is constrained to have no transverse momentum in the lab frame (Fig. F7p). The 



dm 



max 



are 0.385 and 0.998 which are both 



gradients — j^ — measured either side of the kink in Fig. 
very close to their predicted values. 

Also, as expected, if the parent can have very large transverse momentum (SPT), the kink is 
even more pronounced (Fig. [7p). In this case the measured gradients below and above the kink are 
found to be respectively 0.214 and 4.878. Both gradients have changed from the values measured in 
the HERWIG events, coming to within about 5% of the limiting values (attainable only with infinite 
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FIG. 7: (CASE 3v) As for Fig. [6] but for the three-daughter decay, and with mo = 500, rni = 100, and 

TO 2 ,3 = 0. 

momenta) predicted earlier. 

C. Single production with three-daughter cascade decay (CASE 3s) 



The gradients on either side of x = m i f° r three-daughter cascade decays (CASE 3s) are similar 
to CASE 3v (direct decays), except that the mass-shell requirement of the intermediate particle 
places an extra constraint, narrowing the range of allowed values of ttim- Provided that none of the 
decays are near-threshold, the range of values of m^ will be significant, and we therefore expect 
to see kinks in CASE 3s for both SPT and ZPT. We show the ZPT, SPT (medium pt) and SPT 
(large pt) cases in Fig. [si for the particular values of ijiq = 500, mi = 100, m,2 = m^ = 0, m^ = 300. 
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FIG. 8: (CASE 3s) These figure show the TOt(x) distributions obtained from phase-space Monte-Carlo- 
generated events containing a single parent (m = 500) undergoing a sequence of two-body decays via an 
on-shell intermediate particle (m^ — 300) resulting in a final-state containing two massless visible daughters 
(to2 = 7713 = 0) and an invisible massive daughter (ra\ — 100). (a) ZPT (zero pt), (b) SPT (medium pt), 
(c) SPT (large pt), 

Neither of the decays is near threshold and so, as expected, the kink is visible even at ZPT. 
D. Pair production with two-daughter decays (CASE 4) 



We conjectured before that the only possible source for a kink in pairs of two-daughter decays 
(CASE 4) is from large boosts of the sleptons resulting in significant differences in (5 cos 9. In Fig. 19b, 
we find that any kink in m™| x is hardly discernible in pair-decays of sleptons in our HERWIG events. 
This is to be expected, since to get events close to the SPT maximum for pair decays, we need 
both parents to be produced with large energy, and this is unlikely to occur at the LHC. 

When the system of the pair of decaying particles is produced at rest in the transverse plane 
(ZPT), there is no visible kink whatsoever (Fig. K)p). If finally we construct an unphysically large 
Pt distribution for the parents, we confirm that a kink could "technically" become visible (Fig. [9H). 
However, we note that even here the kink is not visible at the yellow-green boundary, and only 
becomes discernable at the green-cyan boundary, indicating that the size of the tail in the m,T2 
distribution which generates the kink is very small - about 1% of the height of the rest of the 
distribution. It seem highly unlikely, therefore, that a kink would ever be seen in CASE 4 in 
presence of any background, not even if the signal cross section was truly gigantic, for detector 
resolution would almost almost certainly smear away the details of this tiny tail. 
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FIG. 9: (CASE 4) As for Fig. [6] but for pair decays, and so plotting the variable rriT^ix)- m o.5 = 300, 
vnifi — 100, m,2j — 0. The outlying points in parts (c) and (d) originate from occasional failed numerical 
minimisations when calculating m^ 2 and have no physical significance. 

E. Pair production with point-like three-daughter decays (CASE 6v) 



It was in a pair of three-daughter decays that the presence of a kink was first noticed pp. 
For our simulations of this case we use the same events as in Section IIVBI but now we treat the 
pair of decays together, recognising that if the parents are pair-produced, then the transverse 
momenta of the invisible daughters are not individually known. For the (HERWIG) SPT case, the 



m.T2 distribution (Fig. 10 1) shows a large number of events near the upper edge. The kink at 
X = rn\ is clear (Fig. |l0p) and the gradients below and above the edge are found to be 0.38 and 
1.00, respectively, to be compared to the SPT global maxima of 0.2 and 5, respectively, and to the 
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FIG. 10: (CASE 6v) As for Fig. [7]but plotting the variable tot2 for a pair of three-daughter decays. mo,i 

500, 7711,6 = 100, and TO2,3,7,8 = 0. 



ZPT maxima of 5/13 ~ 0.38 and 1, respectively. Note that, despite the presence of non-zero pt, 
we find that the observed values correspond very closely to the ZPT maxima, suggesting that the 
relevant events are very close to threshold. 

For the ZPT case we also generate events each containing a pair of virtual three-body decays, 
but this time we use the toy Monte Carlo to impose the requirement that the sum of the parent 



particles has zero transverse momentum Fig. 10 ;. Again we see a kink. 

Notice that the CASE 6v results contrast with those from CASE 4. In CASE 4 a kink was only 
seen in SPT, whereas in CASE 6v kinks are seen in both SPT and ZPT. This means that events 
containing pairs of three-body decays stand a better chance of generating observable kinks at the 
LHC than do events containing pairs of two-body decays, for the former can generate kinks without 
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FIG. 11: (CASE 6v Mtgcii) As for Fig. 10 but for the kinematic variable Mtgcu, rather than witi- 
toq 5 = 500, mi 6 = 100 and m-i 3 7 8 = 0. 



needing to gain large transverse boosts from ISR, while the latter need these boosts. 



Note that the generation of Fig. 10 (and Fig. 12 ) required the Monte Carlo truth record to be 
inspected in order that the quarks could be combined together only in the correct pairings. In 
reality, the pedigree of each jet will be a mystery, and so it is reasonable to wonder whether the 



kink structure of Fig. 10 will remain in a plot where we are blinded with respect to the "truth". 
In m it was claimed that in at least two examples a kink was seen survive the imposition of a 
set of cuts designed to place jets heuristically into the appropriate pairings. Here, however, we 
note an alternative approach which may be followed. Rather than attempt to assign quarks to 
the appropriate parings, one can instead use the MxGen variable [15j in place of rriT2- The MrGen 
variable does not need to be told which momenta come from which side of the event - it treats all 
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FIG. 12: (CASE 6s) As for Fig. [8] but plotting m*T2 for pairs of decays, too, 5 = 500,7711,6 = 100,7712,3,7,8 
0,7774,9 = 300. 



visible momenta on an equal footing. Internally, MxGen performs a minimisation over a (possibly 
large) number of rriT2 evaluations, with an evaluation for each distinct partition of the visible 
particles into two sets (one for each parent). It was shown in [15J that the endpoint of M^Gen 



shares the same properties as the endpoint of rriT2 at X = rn i- I n Fig- H therefore, we look for a 
kink in the maximal values of Mxcen as a function of \t an d indeed the kink is found to be present 
at both ZPT and SPT. Moreover, it is interesting to note that variation with x °f the -^TGen 



maximum seen in Fig. 11 is surprisingly similar to the variation with x of the rriT2 maximum seen 



in Fig. [TOj This is remarkable as it suggests that the endpoint of the M^Gen variable appears to be 
able to capture almost the same information as the endpoint of the txiti variable, even though the 
former is always blind to the correct quark assignments, and the latter expects to be given them. 
Further study is needed to determine which of these two variables is the better to pursue. 



F. Pair production with three-daughter cascade decays (CASE 6s) 



One may conduct the same distributions as in Section IV E but for pairs of cascade decays into 
three particles (CASE 6s) rather than for pairs of point-like three-body decays. The results for 



CASE 6s (Fig. 12 ) are qualitatively the same as those for CASE 6v: a kink is seen in each of the 



SPT and the ZPT distributions. As before, the kink becomes more pronounced as the amount of 
Pt increases: ZPT — > SPT (medium pt) — > SPT (large pt). 
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V. CONCLUSIONS 

We have examined a wide range of cases in which heavy particles are produced singly or in 
pairs at hadron colliders and decay to heavy invisible daughters (as exemplified by the diagrams in 
Fig. FT] and Fig. pi). We have examined how the endpoints of suitable kinematic variables depend on 
the a priori unknown mass of the invisible daughter, paying particular attention to the behaviour 
around its true mass. 

We found that, unless the mass of the invisible particle is known exactly, the global maxima (i.e 
the true endpoints) for the general case (SPT) always correspond to events in which the transverse 
momenta of the particles are asymptotically-large compared to their masses. For x < m ii the 
global maximum occurs when the invisible daughters are emitted in the same direction as the 
velocity vectors of their parents in the lab frame. For \ > m i> t ne global maximum occurs when 
the invisible daughters are emitted backwards relative to their parents' lab velocities. We have 
shown that, because asymptotically large energies are required, it is rather unlikely that significant 
numbers of maximal (or rather, near-maximal) events will be obtained in experiments at energy 
scales not far above the masses of the particles involved. But we have also demonstrated that 
non-maximal events still generate a kink. 

We found that the maximum is different if the extra ZPT condition is applied. The ZPT global 
maximum for \ > m i occurs in events in which the invariant mass, rajy, of the visible iV-daughter 
system takes its maximum value. The global maximum for \ < m i occurs when m^ takes its 
minimum value. For CASE 2 , the invariant mass is fixed, and so in the absence of ISR (i.e. in the 
ZPT case), this is the one case in which the kink does not appear even in principle. In all other 
single-parent decays, the global maxima generate a kink both at SPT and ZPT. 

We note that the ZPT kink has a different origin to the SPT kink. The former is related to 
a degree of freedom possessed by the invariant mass(es) of the set(s) of visible daughters of each 
decaying particle. The latter is an artifact caused by invariant masses being evaluated in a frame 
which is not the parent's rest frame when the invisible particle's mass is hypothesized at incorrect 
values. This second effect is only significant when the parent's velocity is comparable to the speed 
of light. Now, it seems likely that the ISR in events of this type at the LHC is likely to be small, but 
non- vanishing [16]. But since the maxima obtained with vanishing ISR are not maxima (not even 
local maxima) with non- vanishing ISR, it is clear that the behaviour with small, but non-vanishing 
ISR is important. 

To assess the implications of this, we considered the values of tut that result from a finite 
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Case 



2 
3v 

3s 
4 

Gv 
Gs 



Some PT (SPT) 



Kinky 
Kinky 
Kinky 
Kinky 
Kinky 
Kinky 



Zero PT (ZPT) Figures 



Smooth 
Kinky 
Kinky 

Smooth 
Kinky 
Kinky 



10 



11 



12 



TABLE I: Summary of whether or not we expect (and observe) kinks in the different cases, (as shown in Fig.n] 
and Fig. pi). "Kinky" indicates cases in which we observe a discontinuity in the gradient of the maximum 
value of an appropriate kinematic transverse mass variable at x = m i i whereas "Smooth" indicates the lack 
of such a discontinuity. 



distribution of events in j3 and cos #, as defined in Section III B In this formalism, ZPT events 
have (3 = 0. As well as reproducing our previous results, and showing that the global maximum of 
the SPT case is unlikely to be reached, we showed that a kink still arises at (mo,wi), even if we 
do not reach the global maximum. In order to generate a kink, one simply needs events that have 
small relative rapidity, and are kinematically distinct either in their values of m^, or in the values 
of {3 cos 9. So a kink will always be generated in a single-parent decay, unless all events have the 
same value of the invariant mass and the same value of f3 cos 9. Similar conclusions are reached 
in the case of identical pair decays. A full analysis of pair decays, including non-identical decays, 
remains to be done. 

If the bulk of events are SPT (are produced with enough transverse recoil), then kinks are 
ubiquitous, occurring in every case shown in Fig. [2] and in any topology that can be represented 
as Fig. Hi a) or Fig. [lib). By contrast, if the bulk of events are ZPT (produced with little or 
no transverse momentum), we found that for CASE 2 and CASE 4, where the invariant mass is 
fixed, the maximal and minimal gradients are the same, and no kink is found. These results are 
summarized for our exemplar cases in Table |Tj 

We conclude by making some remarks about the experimental viability of the method. Firstly 
we recognise that backgrounds, both from Standard Model and competing new-physics processes 
will need to be considered. Secondly we point out that, in real physics events, the finite widths of 
the unstable particles will have the effect of smearing the end-points. 

When backgrounds, detector smearing and finite width effects are included, it will clearly become 
impossible to determine the position of the endpoint from individual maximal events. Unless the 
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shape of the distribution near the endpoint is known then one must try to infer its position by 
looking for a region of rapidly decreasing numbers of events. The presence of a significant fraction 
of events near the kinematic limit therefore becomes much more important when such effects are 
considered. Since we expect that cases in which kinks only exist only at SPT (i.e. CASE 2 and CASE 
4) will have a small fraction of events potentially contributing to the kink, we have less confidence 
in these cases being experimentally accessible with this method compared to cases which also have 
kinks at ZPT (i.e. all the other cases). 

One might have reasonable confidence that, even after the inclusion of backgrounds, measure- 
ment imprecision, and finite widths, the method presented would still be very useful. After all, 
as remarked in Section [ITJ the Tevatron VF-mass measurements demonstrate that, with adequate 
Monte Carlo modelling of these complicating physical effects, endpoint techniques can be extremely 
useful at extracting mass information. 
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